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We discuss specific heat and neutrino emissivity due to direct Urea processes for quark mat- 
ter in the color superconductive Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) phase of Quantum- 
Chromodynamics. We assume that the three light quarks u, d, s are in a color and electrically 
neutral state and interact by a four fermion Nambu-Jona Lasinio coupling. We study a LOFF state 

' characterized by a single plane wave for each pairing. From the evaluation of neutrino emissivity 

and fermionic specific heat, the cooling rate of simplified models of compact stars with a quark core 

, in the LOFF state is estimated. 
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I. INTRODUCTION 



Neutrino emission due to direct Urea processes, when kinematically allowed, is the most efficient cooling mechanism 
for a neutron star in the early stage of its lifetime . After a very short epoch, when the temperature of the compact 
star is of the order of ~ 10 11 K and neutrinos are trapped in the stellar core @], the temperature drops and neutrinos 
. are able to escape. However, even for smaller temperatures, e.g. below 10 9 K, the direct nuclear Urea processes 
' n — > p + e + D e and + p — > n + v e , which would produce rapid cooling, are not kinematically allowed, because 
energy and momentum cannot be simultaneously conserved. Therefore only modified Urea processes can take place, 
where a bystander particle, present in the reaction, allows energy-momentum conservation. The resulting cooling is 
less rapid because neutrino emission rates turn out to be e v ~ T 8 , much smaller than the emission rate e v ~ T 6 due 
C| ■ to direct Urea processes. 

ij These considerations apply to stars containing only nuclear matter, ff hadronic densities in the core of neutron 
stars are sufficiently large, the central region of the star should consist of deconfined quark matter [|| (we do not 
consider the case of pure quark stars, see e.g. Therefore direct Urea processes involving quarks, i.e. the processes 
d — > u + e~ + v e and u + e~ — > d + v e , may take place and largely contribute to the cooling of the star, ft has been 
shown by fwamoto |5j] that quark direct Urea processes are kinematically allowed and the corresponding emission rate 
for massless quarks is of the order a s T , where a s is the strong coupling constant. This result assumes that quark 
matter is a normal Fermi liquid. However, since the temperature of aged compact stars is sufficiently low, deconfined 
quarks in the stellar core are likely to form Cooper pairs and quark matter could be in one of the possible Color 
Superconductive (CS) phases (see [||, Q and for reviews This is due to the fact that the critical temperature of 
CS matter is of the order of dozens MeV, well above the estimated temperature of the stellar core < f 00 KeV. 

At asymptotically high densities the energetically favored CS phase is the color-flavor locked (CFL) phase, in 
which light quarks of any color form Cooper pairs with zero total momentum 0, and all fermionic excitations are 
gapped. The corresponding neutrino emissivity and specific heat C are suppressed by a factor e _A / T , where A is 
the quasiparticle gap in the CFL phase. Therefore the cooling of quark matter in the CFL phase is distinctly less 
rapid than in the normal phase. In the CFL phase quark masses can be neglected and color and electric neutrality 
conditions are automatically implemented. However at densities relevant for compact stars the quark number chemical 
potential y, cannot be much larger than 500 MeV and effects due to the strange quark mass m s must be included. 
Requiring that bulk quark matter is in weak equilibrium and electrically and color neutral 0, llfl fllL Ha . , together 
with m s 7^ 0, determines a mismatch Sfi between the Fermi momenta of different quarks, with Sfi depending on the 
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in- medium value of m s . For values of m s less than a critical value, CFL is the energetically favored phase. For larger 
values the CFL phase cannot be realized and quark matter should pair with a less symmetric pattern. 

The ground state of quark matter in these conditions is still a matter of debate and several possible superconductive 
phases have been suggested (see for a review). Recently two superconductive phases characterized b y s apless 
fermionic excitations, i.e. the gapless-2SC (g2SC) phase lip an d the gapless-CFL (gCFL) phase \wi UJl 1is[ have 
been largely discussed. However, it has been shown 0, 12(1 l2ll 12^ that both are "chromo-magnetically unstable" 
because the Meissner masses of some of the gluons associated with broken gauge symmetries are imaginary. 

Another pos sibility that has attracted theoretical attention is the Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) state 
|23l I24L I2H |26|. where the total momentum of the pair does not vanish and counter propagating color currents are 
spontaneously generated. A simplified ansatz "crystal" structure with 

3 

(ip ai (x)C'y 5 'ip/3j(x)) oc ^2 A/ e 2,qrr e ttj 8jey.r , (1) 
1=1 

— 1,2,3 flavor indices, a, /3 = 1,2,3 color indices) has been studied in Refs. [27l |28| and found energetically 
favored with respect to the gCFL and the unpaired phases in a certain range of values of Sfx. In Eq. (JIJ 2 qi represents 
the momentum of the Cooper pair and the gap parameters Ai, A2, A3 describe respectively d — s, u — s and u — d 
pairing. For sufficiently large [i the energetically favored phase is characterized by Ai = 0, A 2 = A 3 and q 2 = q3. 
Gluon Meissner masses corresponding to broken color generators have been evaluated in [2!j and thisphase results 
to be chromo-magnetically stable. More complex crystal structures have been recently proposed in [3(|. From the 
evaluation of the corresponding free-energy in the Ginzburg-Landau approximation one finds that crystal structures 
with more plane waves are energetically favored with respect to the normal phase and the gCFL phase in a wider 
range of densities. We note that by (1) we assume attractiveness in the color antisymmetric channel. This follows 
from the one gluon exchange diagram of QCD. It dominates the asymptotic regime and we assume that it favors 
attraction also at moderate densities (/1 ~ 500 MeV). In other words, though at moderate density nonperturbative 
effects can play a role we are assuming that the superconductive ground state is qualitatively similar to that of very 
high density. 

The aim of the present paper is to evaluate the neutrino emission rate and the specific heat of quark matter in the 
LOFF superconductive phase. We show that, due to the existence of gapless modes in the LOFF phase, a neutron 
star with a quark LOFF core cools faster than a star made by nuclear matter only. This follows from the fact that in 
the LOFF phase neutrino emissivity and quark specific heat are parametrically similar to the case of unpaired quark 
matter (e v ~ T 6 and C ~ T respectively). Therefore the cooling is similar to that of a star comprising unpaired 
quark matter. Incidentally we note that the dominance of quasiparticle gapless modes has been demonstrated also 
for the g2SC phase [3j| and for the gCFL phase [32] , though in the latter phase the behavior of neutrino emissivity 
and quark specific heat is different [e v ~ T 55 and C ~ T 5 ). For calculation of neutrino emissivity in other models 
see [33j . 

Our results should be considered as preliminary, since, as we have already pointed out, the simple ansatz Q should 
be substituted by a more complex behavior as in |30| . If the condensate is the sum of more plane waves, the calculation 
of e u and C is more involved, basically because one should obtain information on the quasiparticle dispersion law 
in the framework of the Ginzburg Landau expansion. For complex condensate patterns resulting from several plane 
waves this is a complicated task (a calculation for two flavors is in j34[) and we leave it as a future work. Given 
these limitations it would be fruitless to consider sophisticated star models. Therefore we evaluate the cooling rate 
by employing toy models of stellar objects, i.e. stars made of nuclear matter with a core in the color superconductive 
LOFF phase. Also, we will not consider the contribution to the cooling of the compact star due to other processes 
such as the neutrino pair bremsstrahlung from nuclei in the crust and pionic reactions |3^] that should also be included 
in a realistic description of the compact star cooling Q . 

Our paper is organized as follows. In Section[n]we briefly discuss the LOFF phase characterized by the condensate 
P. In Section ITTT1 we evaluate the neutrino emissivity and in Section IIVI the specific heat in the LOFF phase. In 
Section[V]we evaluate the effect of l//z corrections, taking into account the results of [3(j . In Section lvTl we estimate the 
cooling rate of toy models of neutron stars. Given our approximations any comparison with experimental observation, 
similar e.g. to those already appeared in the literature [U IH, |3^, H3> HE IS H^, is premature and we limit our 
analysis to a qualitative comparison among various simple models. In Section [VIII we draw our conclusions and in the 
Appendix we report the dispersion laws that are used to compute the quark mixing coefficients and gapless points in 
the LOFF phase. 
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II. NEUTRAL LOFF QUARK MATTER 

Non interacting quark matter consisting of massless u and d quarks and s quarks with an in-medium mass m s can 
be described by the Lagrangian density 

Co = & Q (i9tf- M tf + tif to) V« , (2) 

where i, j — 1,2,3 are flavor indices and a, /3 = 1,2,3 are color indices; Dirac indices have been suppressed; the mass 
matrix is given by M"P = 5 al3 diag(0, 0, m s )ij and = dS al3 Sij. The quark chemical potential matrix is as follows: 



with Qij = diag(2/3, — 1/3, — l/3)y the quark electric-charge matrix and T 3 and T g the Gell-Mann matrices in color 
space. We are interested in the region close to the second order phase transition point between the normal phase and 
the Color Superconductive phase, where, to the leading order approximation in Sfi/n, /i3 = Ms = and jj, e = m 2 s /A^, 
as in the unpaired phase [3(| . In all the numerical estimates we use the value [i — 500 MeV. 

In order to describe quark interaction we employ a Nambu-Jona Lasinio (NJL) model in a mean field approximation. 
The use of the NJL model can be motivated by renormalization group analyses [44J showing the dominance of the four- 
fermion interactions and in particular the numerical importance of the coupling mimicking one gluon exchange. We 
assume here that these results hold both at high and moderate hadronic densities. The resulting quark condensate 
is antisymmetric in color and flavor indices and we assume the behavior (Q for it. This phase has been studied 
employing a Ginzburg-Landau (GL) expansion in Ref. l27l . Requiring color and electric neutrality results in Ai = 
and, to the leading order in 1/fi expansion, A2 = A3. In |3|j it has been shown that such corrections affect this result 
so that A2 < A3. These effects will be discussed in Section IVI while in this Section we take into account only the 
leading order effect in Sfi/fj,. 

The GL approximation is reliable in a region close to the second order phase transition point [2^] where the favored 
"crystal" structure is characterized by q2=q3=q. The LOFF phase is energetically favored with respect to the gCFL 
phase and to the normal phase in the range of the chemical potential mismatch y £ (130, 150) MeV, with 



Moreover the magnitude of q is given by 



(4) 



(5) 
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and A = A 2 = A 3 < 0.3A where A is the value of the gap in the two-flavor homogeneous case [3(j. We use the 
value Ao = 25 MeV. Recent analysis [3^ have shown that the range of values of y where the LOFF state prevails 
can be enlarged assuming a more complex pattern of space dependence for Aj (r) . For the purposes of this paper it 
is however sufficient to consider the simple ansatz ijljl. 

In the basis A = (1, 9) = (u r , d gi Sb, d r , u g , s r , Ub, s g , db) the quark chemical potentials can be expressed by 

fi A = M + MA (6) 

with 

y V - % _ y _y ,„n 
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Here 



z = cost? , 



(9) 
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with $ the angle between the quark momentum and the pair momentum. For later convenience we also define 

/i A (0) =n A (q = 0) . (10) 

We note explicitly that the previous chemical potentials can be obtained by a redefinition of the quark fields in the 
following way: 

4> ai (x) = e l ^ x ^ m (x) , (11) 

where q a i = q for (a, i) = [(r, d), (g, u), (r, s), (6, u)], q„ = 2q for (a, i) — (r, u) and q„ = in the other cases. Such 
a redefinition of the quark fields allows to eliminate the space dependence of the condensate. However quark momenta 
are shifted as follows: 

Pn r ^P« r -2q, Pd g ->Pd g , P s „ -> P Sb , 

Pd r Pd r - q , Pu g -> Pu g - q , p Sr -> p Sr - q , p« b -> p« b - q , 

Ps 3 -> Ps 3 , Pd b -> Pd b • (12) 

III. NEUTRINO EMISSIVITY 

The transition rate for the /3 decay of a down quark d a , of color a = r,g, b, into an up quark u a 

d a (pi) -> ^efe) + w Q (p 3 ) + e _ (p 4 ) (13) 

is 



Wfi = V(2n) i S%p 1 -p 2 -ps-p i )\M\ 2 l[— , (14) 



i=l 



where V is the available volume and M. is the invariant amplitude. Neglecting quark masses the squared invariant 
amplitude averaged over the initial spins and summed over spins in the final state is 

\Mf = 6AG 2 F cos 2 9 c ( Pl .^(ps • p 4 ) , (15) 

where Gf is the Fermi constant and 9 C the Cabibbo angle; we will neglect the strange-quark f3 decay whose contribution 
is smaller by a factor of tan 2 6 C in comparison with (|15l) . For relatively aging stars, there is no neutrino trapping 
which means that neutrino momentum and energy are small. The magnitude of the other momenta is of the order of 
the corresponding Fermi momenta p F , p F and pp. As discussed in Sectional f° r d and u quarks they are of the order 
of 500 MeV; the electron momentum is of the order /j, e , which is smaller, but, due to the assumptions discussed in the 
previous section, still sizeable. On the other hand the neutrino momentum p 2 is of the order ksT. It follows that the 
momentum conservation can be implemented neglecting p2 and one can depict the 3-momentum conservation for the 
decay (|13|l as a triangle 5] having for sides pi, P3 and p 4 . It follows that we can approximate 

(Pi -P2)(P3 -Pi) - E 1 E 2 E 3 E i (l - cos 0i 2 ) (1 - cos# 34 ) (16) 

where Ej are the energies and #12 (resp. #34) is the angle between momenta of the down quark and the neutrino (resp. 
between the up quark and the electron). Neutrino emissivity £„ is defined as the energy loss by neutrino emission 
per second per unit volume. To compute it, we have to multiply (|15|l by the neutrino energy and thermal factors, 
integrate over the available phase space and sum over the three colors 0. Moreover, since the states having definite 
energy are the quasiquarks involved in the Cooper pairs, one has to multiply the result by two mixing (Bogoliubov) 
coefficients depending on the quasiparticle dispersion laws 32]. Therefore we get 



a=r,g,b a—r,g,b 



n 



d z p l 



E 2 W fi 



(2tt) 3 

■n( Pl ) • [1 - n( P3 )] ■ [1 - n( P4 )] • B 2 da ( Pl )B 2 Ua (p») (17) 
where n(pj) are thermal equilibrium Fermi distributions: 

^ TT^r ■ (18) 
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*, = • (19) 

They are appropriate here because strong and electromagnetic processes establishing thermal equilibrium are much 
faster than weak interactions. The overall factor of 2 keeps into account the electron capture process. 

Thus far the results are similar to those of jj^]. However, in the present case, the phase space integrations and the 
momentum conservation law must be treated with great care. As a matter of fact, in the LOFF color superconductive 
state d a and u a are in general paired with quarks of another color, and there is breaking of translational invariance 
because the total momentum of the Cooper pair is 2q. Let us choose the z— axis along the direction of q and let dj 
and 4>j be the polar angles of pj. Since the temperature is much smaller than the gap parameter, the dominant modes 
in the d— and momentum integrations are the gapless ones, that we denote as Pj (with j = 1 for down quark and 
j = 3 for up quark). They arc defined as follows: Pi (resp. P3) is the quark down (resp. quark up) momentum where 
the corresponding quasi-particle energy vanishes (for more details see the Appendix). On the other hand the relevant 
momentum for the electron is its Fermi momentum. Therefore we have 

d 3 pi J d 3 p 3 J d 3 p A fa J fildpidQi P? dp x d^P 2 dp 3 dfl 3 , (20) 

with dflj = sin $j ddj dipj . While these approximations are similar to the assumptions used in [S2l | in dealing with 
the homogeneous gCFL phase, in the LOFF phase the gapless momenta Pi and P3 depend on the angle $j that 
quark momenta form with the pair momentum 2 q. In order to simplify the expression of the integral in Eq. I|17|) the 
variables pj can be traded for Xj, as in Eq. (|19J) . We put X2 — E2/T for the neutrino, X4 = (p — /i e )/T for the electron 
and Ej = [ij + £j(p) for the quarks. Expanding around the gapless modes one has Ej{p) ~ fij + Vj (p — Pj) (for j = 1 
and 3), with the quasiparticle velocity given by 



_ 

3 dp 



(21) 



As discussed in detail in the Appendix the dispersion law of each quasiparticles has from one to three gapless modes. 
Therefore one has to expand the corresponding dispersion laws around each gapless momentum. Another point to 
be stressed concerns the conservation of three-momentum. As we have noted, to get rid of the space dependence 
in the condensate, one redefines the quark fields, see Eq. which amounts to a redefinition of momenta as in 

Eq. JI3). Using the new quark fields in the weak decay matrix element (115(1 adds a momentum — q to the momentum 
conservation law that now reads: Pi — P3 — P4 — q = (where the neutrino momentum has been neglected). 

Employing the above approximations the neutrino emissivity for each pair of gapless momenta Pi , P3 , can be written 

as 

327r s J (27r) J ^J- J \vi \ \v 3 \ 



G 2 F cos 2 9 c ^ e T\ f dv j-j PlPiBj^BljP,) 

3 = 1,3,4' 



Stt* J (2tt)3 11 / > MH 



x U Pj+^-^Pl^l-Pj-^ e - 4r .( Pl - P 3-p 4 - q ) 7 (22) 

where P\ and P3 depend, respectively, on the angles $1 and $3 and E>i a is a mixing coefficient arising from the 
Bogoliubov-Valatin transformation representing the probability amplitude that the gapless quasiparticle has flavor i 
and color a. Moreover Hi liFfl 



/+00 /-+oo r+<x> n+00 45 77T 6 

dxi / dx 2 x\ / dx 3 / dx4n(xi)n(-x 3 )n(-X4)5(xi - x 2 - x 3 - X4) = — — . (23) 

00 Jo J-00 J-00 5040 

Some of the angular integrations appearing in (|22|l can be performed analytically, with the result 



J />+OC />1 

lY\\ drr 2 d{cos-d)e lqrmai> 
k=1 Jo J-i 



G 2 F cos 2 e c ji\ T 6 
4ir 6 

x / d(cos^i)e- tPircos ' ?cos1?l J (Pirsini?sin-(?i)/ fc (Pi) 1 11 
Juj n \vi\ 
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"3 p2 q1 lp\ 

d{cosd 3 )e tP3rcos ^ cosi>3 J (P 3 rsmdsm^ 3 )g k {P 3 )- 3 " ' 



^2 
1 



x / d(cos^ 4 )e v ' !rcos ' 5cosl ' 4 Jo(/i e rsin^sin^4) 



P\ + q 2 , fs = - 2Piqcost?i, and g x = 1, g 2 = 55 = - 



(24) 

Here for convenience we have decomposed the integral as a sum over various terms, with fi = f 3 = fi = X, fi = 

2 Me P 3 ' 53 2//, ' ff4 2P 3 ' 
In Eq. I|24[l J is the Bessel function of zeroth order; u>q, U\, u>2i w 3 are appropriate limits taking into account 
kinematic constraints; "& is the angle between the directions of r and q. Even if each quasiparticle dispersion law 
is characterized by various gapless momenta, the number of the gapless momenta relevant for the evaluation of the 
neutrino emissivity can be reduced observing that the momenta pi, P3 and P4, with |p4| ~ /i e , must satisfy the 
condition pi — P3 — P4 — q = 0. 

We now consider the contributions of the three colors separately. We discuss in detail only the decay of the blue 
quark that also gives the largest contribution. The decay of the other colors is treated analogously. The blue down 
quark is unpaired because in the Ginzburg-Landau approximation the gap parameter Ai vanishes |2?| . Therefore, 
as discussed in the Appendix, the dispersion law is Ei(p) = p, with gapless momentum P\ = mixing coefficient 
S d6 (P 1 ) = landi; 1 = l. 

Let us now consider the up-blue quark. It is paired with the strange-red quark with gap parameter A 2 = A. In 
the Appendix fsubsection lA 3|) we give the dispersion law and the gapless momentum: 



Pi 




M6 + 



A 2 . 



(25) 



We do not consider the momentum corresponding to the other gapless mode of Eq. I|A21|I . as it does not satisfy the 
condition of momentum conservation. The emissivity for the decay of the blue quark is therefore: 

r-l 



where 



_bluc 



G\ cos* c £i4l* l [+°° 



d(cos$i)e 



—iP\r cos cos $1 



d(cosl?)e i9rt 



Jo (Pi?" sin $ sin $1) 



d(cos i? 3 ) e lPar cos * cos 1,3 J (P 3 r sin d sin tf a )g (cos i? 3 , cos #i ) 
d(cos #4 ) e ^ r cos * cos 1,4 J (fJ. e r sin 1? sin 4 ) , 



LU 2 



Mm < z, 



8A 

y 



(26) 



(27) 



and 



g(cosd 3 , cost?i) = ^1 — 

x e [ 1 



,i 2 + q 2 - 2qfx d cos 1?! - P 2 - ,i 2 \ [P 3 P Mb (P 3 )] 



2P 3 /Xe 

2(7/1^ cosi?! 



P 2 
^3 



«3 



«3 



P.3- 



2P 3 fi. 
Me + M7 



P.3- 



^6 + M7 



A 2 



(28) 



(29) 



We note that the presence of L02 implements the existence of the gapless momentum P 3 while the Heaviside function 
implements the condition | cos $3 1 < 1 . 

The case of quarks with colors red and green can be treated in a similar way, though the numerical computation is 
more involved because for these colors neither the down nor the up quarks are unpaired. Because of this one expects 
that the emissivity of these quarks is smaller than for blue quarks, an expectation confirmed by the numerical analysis. 
We will report our numerical results in Section I VII 
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IV. SPECIFIC HEATS 

At low temperatures the largest contribution to specific heat C is determined by the sum of the specific heats of 
the fermionic quasi-particles. In the three flavor LOFF phase the quasi-particle specific heats are given by 

^ ej^, (30) 



(2tt) 3 j dT 

where j = u r ,d g , Sb ,d r ,u g , s r ,Ub , s g , db and 6j = \Ej (p) — [ij | are the quasi-particle dispersion laws. Since we work 
in the regime T <C A <C //, the contributions of gapped modes are exponentially suppressed and each gapless mode 
contributes by a factor oc T. This results follows from the evaluation of the integral in Ea. (|30|l employing the saddle 
point method and assuming that the quasi-particle dispersion laws are linear in the gapless momenta. For the present 
choice of parameters, the quasiparticle dispersion law is quadratic in the gapless momenta in a negligible range of 
values of the angle between the direction of quasiparticle momentum and the direction of q. Such a situation is 
quite different from the homogeneous gCFL case where a quasiparticle dispersion law is quadratic 32] and gives the 
dominant contribution to the specific heat. 

Within the above-mentioned approximations the fermionic specific heath is given by 



C ^3 



J 2 l\v I \ + \v n \ + \vm\ + + K bSr \ \&*\ ( } [ ] { ] 



where the various gapless momenta and Fermi velocities on the r.h.s have been defined in the Appendix and in Eq. 
(|21|l . In particular Pj, Pjj and Pm are the gapless momenta of the sector (u r ,d g , Sb) and Vi, vn and vjjj the 
corresponding velocities. We have also added the electron contribution, though /i e is numerically much smaller than 
all the quark momenta P g (generally of the order of We also note that the integration range is in general smaller 
than (—1, +1) because the quasi-particle dispersion laws are gapless only in a restricted range of values of the variable 
cos??. 

V. ROLE OF 1/n CORRECTIONS 

In previous sections we have evaluated emissivity and specific heat in the three flavor LOFF phase of QCD. All 
the calculations are based on the High Density Effective Theory, see ^B,E3;E3|, where one takes the limit fj, — * oo. 
In such approximation one neglects the contribution of the antiparticles; moreover the m s effects are treated at 
the leading order by a shift in the strange Fermi momentum pf w /i s — m 2 s /2fi. One can show |36| that, as already 
stressed in Section [H] in this approximation the electron chemical potential /i e is given by [i e = m^/A/j,, which implies 
a symmetric splitting of the s and d Fermi surfaces around the u Fermi surface. Therefore Ai = and A2 = A3. If 
\x is not large enough this approximation is not justified and higher order corrections must be included 

The next-to- leading corrections were computed in [36|. One expands the strange quark momentum up to the 
next-leading order: 

Then one substitutes this expression in the Lagrangian in Eq.Q, getting a correction. It can be proven |3r| that 
other corrections, e.g. the antiparticle contribution, is next-to-next-to-leading in the weak coupling approximation 
and therefore can be safely neglected. 

Keeping only terms of 0(A/q) 4 in the free energy, the 1/fi shift in Eq. (|32|l results in the chemical potentials 

2 A 

m mz 

Me= V~4V' ^' M8 = - (33) 

The introduction of the 1 j [i terms in the free energy results in an asymmetric splitting of the Fermi surfaces of the 
u, d and s quarks. As a matter of fact one gets in the normal phase 

Mu " pf _ m 2 5 m\ - fi u = m 2 s 1 m\ 

2 8^ 96 /i 3 ' 2 8^ 96 fi 3 1 ' 

As a consequence one expects in the LOFF phase still Ai = 0, but A2 < A3. For example for fj, = 500 MeV one finds 
at rn 2 /fi = 140 MeV that A2 = and A3 ~ 0.35Aq. Thus for this value of m s the vacuum consists of LOFF pairs 
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of u r — d g and u g — d r quarks and Fermi seas of unpaired s, Ub and quarks. This new phase was called L0FF2s 
in [33 since it is a two flavor LOFF phase whose excitation spectrum consists of quasi-particles with dispersion laws 



M4 - Ms 




"4 + «e 



A 2 
^3 



(35) 



and a sea of unpaired blue and strange quarks. 

Motivated by these results we study neutrino emissivity and specific heat of the LOFF2s phase. Since blue up and 
down quarks are unpaired, their contribution to the neutrino emissivity is 1/3 of the neutral unpaired quark matter. 
Using for it the estimate of Ref. @ one obtaines 



_blu 
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- X 

3 



457 G% cos 2 e r 



630 



a s fJ, u fidfJ-e (k B T) e 



(36) 



where we have restored the correct factor of c, h and ks- We note that in this equation a s = a 1 /4tt = 4a c where g 
gives the coupling among quarks and gluons and a c is the coupling constant appearing in Ref. Jfj . Notice that the use 
of this estimate by Iwamoto @ is only indicative, because at the densities we are considering, perturbative QCD is 
unreliable. Therefore (|36|l can provide at best an order of magnitude estimate. Neglecting the contribution of strange 
quarks, the other contribution to the emissivity is determined by red and green light quarks and can be computed 
following the same lines of Section ITTT1 Since the LOFF pairing results in a restriction of the available phase space, 
the red and green channels contribute to the total decay rate less than the blue one, similarly to the discussion in 

Section urn 

Next we turn to the specific heat, which is given by 



C 



2 x 



0P+) 2 , (P-) 



+ (^) 2 + (^) 2 + {V Ub ? + (/ b ) 2 + {^f + Me 



(37) 



where the gapless momenta P± are defined as 



p±=M 4 + M 5± jM4-M5_ A 2 



(38) 



and the Fermi velocities are given by 



v + = V- 



M4 - M5 



A 2 
^3 



M4 - Ms 



(39) 



The factor 2 in the first two addenda of the r.h.s. of Eq. Ij37(l takes into account the fact that the dispersion laws for 
the u r — d g quasiparticles are the same of the u g — d r ones. 



VI. COOLING BY NEUTRINO EMISSION 



Let us assume the presence of quark matter in the LOFF state in the core of a compact star. This will affect the 
cooling process and we now discuss these effects by comparing various models of stars along lines similar to Ref. [13 . 
Given the approximations used in the study of the LOFF phase, in particular the use of the simple ansatz instead 
of more complex space behavior of the condensate, it would be fruitless to employ sophisticated models. Therefore we 
will use a simplified approach based on the study of four different star toy models. The first model (denoted as 1) is 
a star consisting of noninteracting "nuclear" matter (neutrons, protons and electrons) with mass M = 1.4Mq, radius 
R = 12 km and uniform density n = 1.5 no, where no = 0.16 fm -3 is the nuclear equilibrium density. The nuclear 
matter is assumed to be electrically neutral and in beta equilibrium. The second model (II) is a star containing a 
core of radius Ri = 5 km of neutral unpaired quark matter in weak equilibrium at fi = 500 MeV, with a mantle 
of noninteracting nuclear matter with uniform density n. Assuming a star mass M = 1.4 Mq from the solution of 
the Tolman-Oppenheimer-Volkov equations one gets a star radius R2 = 10 km. Finally we discuss two simplified 
models of compact stars containing a core of electric and color neutral three flavors quark matter in the LOFF phase, 
with jji = 500 MeV and m 2 / fx = 140 MeV and Aq = 25 MeV. Both models have a mantle of noninteracting nuclear 
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matter and differ by the values of the LOFF gaps. In model III we neglect 0(1/ n) corrections and we take Ai = 
and A 2 = A 3 = 0.25 A , the value of the gap parameters in the LOFF phase at y = 140 MeV This is the 

approximation discussed in Section [H] In model IV we include the 1/fi correction discussed in Section Ivl and we take 
the values of the gaps as in the LOFF2s model discussed in Ref. j3||, i.e. Ai = A 2 = 0, A3 = 0.28 Ao for m^/fi = 140 
MeV. Since the values of the gaps in both LOFF models are small the radii of the star and of the quark core do not 
differ appreciably from those of a star with a core of unpaired quark matter, i.e. R\ — 5 and R2 = 10 km (also in 
these cases M = 1.4M Q ). 

The main mechanisms of cooling are by neutrino emission and by photon emission, the latter dominating at later 
ages. Therefore the star cooling is governed by the following differential equation: 



dT _ L v + L 7 _ V nm e™ + V^eg" + 



dt V r nm + V r qm V r nm + V r qm 



(40) 



Here T is the inner temperature at time t, while L v and L 7 are neutrino and photon luminosities, i.e. heat losses 
per unit time. Neutrino luminosity is obtained multiplying the emissivity by the corresponding volume. Therefore 
one must distinguish between the emissivity e" m of nuclear matter, which is present in the volume V nmi from the 
emissivity due to quarks, if they are present in some volume V qm ; Cy m and Cy 1 denote specific heats of the two forms 
of hadronic matter. 

For the emissivity of nuclear matter e™ m we use the standard value 



"0 



2/3 



T 



1.2 x 10 4 erg cm^s" 1 ) ( — ) ( ) (41) 



arising from the analysis of the modified Urea processes n + X—>p + X + e + D, with X — p or n. As for the quark 
contribution e^ m , we consider both unpaired quarks and quarks in the LOFF state, depending on the model. The 
former contribution is denoted £ q m < un P a ' lred and is given by Ea. i|3t)[l multiplied by a factor of 3. The latter have been 
estimated in Sections IIIII and Ivl for the two models of LOFF quark matter. For all these stars we assume a common 
interior temperature T, since the matter comprising compact stars is made of good conductors (p, n or quarks). 

Let us now turn to cooling by photon emission, the dominant process for sufficiently old stars (t > 10 6 years). In 
this case the luminosity can be estimated in the black-body approximation as follows: 

L 7 = 47ri? 2 crT s 4 , (42) 

where R is the radius of the star, a is the Stefan-Boltzmann constant and 

/ \ ^ 4 / T \ °' 55 

T -(°- 87xl ° 6K ) 4 (i()4^) (wk) ' (43) 

is the surface temperature |49l. l5(ij| . with g s = GjyM/R 2 the surface gravity. 

Let us now discuss specific heats. They are given for the LOFF model III and for the LOFF2s model IV by Eqs. i|31|) 
and l|37fl respectively. For nuclear and unpaired quark matter the total specific heat is given by the sum of the specific 
heats of the different fermionic species. They can be approximated by the fermionic ideal gas result 



Cy 



31 

3h 3 c 



PF\/m 2 c 2 + (p F ) 2 , (44) 



where m is the fermionic mass and pf the Fermi momentum. For non-interacting nuclear matter, the three species 
are neutrons, protons and electrons with Fermi momenta evaluated as in neutral matter in weak equilibrium |l|: 

p n F = (340 MeV) ( — 

P P f = Pf = (60 MeV) f — J . (45) 
For neutral unpaired quark matter in weak equilibrium, there are nine quark species, with Fermi momenta independent 

9 1 r 2 

of color and given by p% = fi+ —, p u F = fi - — , p s F = fi - 

Eq. I|40|) is solved imposing a given temperature To at a fixed early time to (we use To — > 00 for to — > 0). In Figs. 
m21and[3]the cooling curves for various models of stars are shown. In Fig. ^the inner temperature, in Kelvin, as 
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FIG. 1: (Color online) Inner temperature, in Kelvin, as a function of time, in years, for three toy models of pulsars. Solid 
black curve refers to model I; dashed line (red online) refers to model II; dotted curve (blue online) refers to model III. Model I 
is a a neutron star formed by nuclear matter with uniform density n = 0.24 fm~ 3 and radius R = 12 Km; model II corresponds 
to a star with R2 = 10 km, having a mantle of nuclear matter and a core of radius Ri = 5 Km of unpaired quark matter, 
interacting via gluon exchange; model III is like model II, but in the core there is quark matter in the LOFF state, see text for 
more details. All stars have M = 1.4 M®. Parameters for the core are fi — 500 MeV and m 2 s /^140 MeV. 



a function of time, in years, for three models is reported. Solid line (black online) is for model I (electrically neutral 
nuclear matter made of non interacting neutrons, protons and electrons in beta equilibrium); dashed curve (red online) 
refers to model II (nuclear matter mantle and a core of unpaired quark matter, interacting via gluon exchange); the 
dotted line (blue online) is for model III (nuclear matter mantle and a core of quark matter in the LOFF state in the 
leading 1/n approximation, as discussed in Section^J. The curve reported in Figs. n an d|3for model III corresponds 
to m s — \/140 (i MeV, however with increasing values of m s the neutrino emissivity decreases. This is due to the fact 
that A decreases as one approaches the second order phase transition to the normal state. On the other hand in this 
case the quark matter in the star tends to become a normal Fermi liquid, for which the description of Iwamoto which 
includes Fermi liquid effects, must be adopted. 

Fig. |21 gives the star surface temperature as a function of time, as obtained by use of Eq. (|43|l . The three curves 
refer to the same toy stars as in Fig. ^ Both diagrams are obtained for the following values of the parameters: 
fx = 500 MeV, ml/fj, — 140 MeV, Ai = 0, A 2 = A 3 ~ 6 MeV. For unpaired quark matter we use a s ~ 1, the value 
corresponding to /1 = 500 MeV and Aqcd = 250 MeV. The use of perturbative QCD at such small momentum scales 
is however questionable. Therefore the results for model II should be considered with some caution and the curve is 
plotted only to allow a comparison with the other models. In particular the similarity between the LOFF curve and 
the unpaired quark curve follows from the fact that the LOFF phases are gapless, so that the scaling laws cy ~ T and 
e v ~ T e are analogous to those of the unpaired quark matter. However the cooling curve for unpaired quark matter 
depends on the value we assumed for the strong coupling constant, a s ~ 1, which is an extrapolation to a regime 
where perturbative QCD is less reliable. Therefore the similarity between the curves of models II and III might be 
accidental. 

In Fig. |21 we compare two models of a star with a nuclear mantle and a quark core in the LOFF state. The 
appoximation used for these LOFF states were discussed in Sections ITU and Ivl The continuous curve (red online) is 
for model III and the dashed black curve is for model IV, i.e. with quarks in the LOFF2s state. One can note that 
both curves for the LOFF models are similar and show a rapid cooling, much faster than for ordinary stars comprising 
only nuclear matter. This implies that the results for models with a LOFF core are rather robust and, at least for 
relatively young stars, the presence of quark matter in the LOFF state should offer a signature clearly distinct from 
that of an ordinary neutron star. 
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FIG. 2: (Color online) Surface te i — -~ m : " T - r - , - : - t— f — : "-^ three toy models of pulsars 

described in Fig. 




Log t [yrs] 



FIG. 3: (Color online) Comparison of the cooling curves for the two toy models of stars with a core of quarks in a LOFF state. 
The curves express the inner temperature, in Kelvin, as a function of time, in years, for model III (solid, red online, curve) and 
model IV (dashed black curve). Model III is a star with a radius R2 = 10 km, having a mantle of nuclear matter and a core 
of radius Ri = 5 Km of quark matter in the LOFF state in the leading approximation (Ai = 0, A2 = A3 ~ 6 MeV); model 
IV is analogous to model III, but the core is in the LOFF2s phase (Ai = A 2 = , A 3 ~ 7 MeV). All stars have M = 1.4M . 
Parameters for the core are fi = 500 MeV and m 2 J ju=140 MeV. 



VII. CONCLUSIONS 



Our study can be summarized as follows. We have computed neutrino emissivity and specific heat for quark matter 
in the LOFF state of QCD in presence of three light quarks. Quark matter has been assumed to be in weak equilibrium 
and in a color and electrically neutral state. We have considered the simplest ansatz for the condensates (single plane 
wave for all LOFF pairings). We have studied two models of the LOFF state, differing in the approximations used 
for the Ginzburg Landau evaluation of the free energy and the gap equation. Our analysis shows a similarity between 
the two approximations, which points to a robustness of our results. We have used this study for an estimate of the 
cooling of compact stars with a nuclear mantle and a quark core. 

Which conclusions can we draw from the present study? From Figs. ^andElwe see that stars with a LOFF core cool 
down faster than ordinary neutron stars. This might have interesting phenomenological consequences. At the present 
time observational results on the cooling of pulsars are being accumulated at an increasing rate (for compilation of 
data and comparison between theoretical models and data sec, e.g. HI S IM HJ EJ 13 Hi). Some data indicate 
that stars with an age in the range 10 3 — 10 4 years have a temperature significantly smaller than what is expected 
on the basis of the modified Urea processes. It is difficult however to infer, from these data, predictions on the star 
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composition, as the stars may have different masses. As for the impact on our study, it is useful to repeat here that 
our analysis should be considered as preliminary because the identification of the quasiparticle dispersion laws for 
the favored crystalline LOFF structures, formed by more plane waves, is still lacking. This is the reason why we 
have not tried a comparison with observations in this paper. Nevertheless some qualitative assessments can be made. 
Quantum Chromodynamics predicts that at the densities that can be reached in the core of compact stars deconfined 
quark matter should be present, and, if so, it should be in a color superconductive state, since Cooper condensation 
of colored diquarks is energetically favored. Slow cooling is typical of stars containing only nuclear matter or of stars 
with a color superconductive quark core in a gapped phase (e.g. CFL). If a careful comparison with the data could 
allow to rule out slow cooling for star masses in the range we have considered, this would favor either the presence 
of condensed mesons or quark matter in a gapless state in the core, since gapped quarks emit neutrinos very slowly. 
Meson condensation also might allow rapid cooling. However reliable calculations based on the chiral effective field 
theory [SJ can only be made for m 2 s ~ /iA, i.e. far away from the region where the LOFF state is stable; therefore no 
direct comparison between the two phases is available. For intermediate densities the quark normal state is less favored 
than color superconductive states for a wide range of values of the baryonic chemical potential fi and the strange 
quark mass m s . This leaves us with gapless quark phases and, among them the LOFF state is favored since we know 
that the gapless phases with homogeneous gap parameters such as, e.g. the gCFL or the g2SC phase are instable, 
while the LOFF phase does not suffer of a similar instability. Let us finally observe that our conclusions should remain 
valid, at least qualitatively, also for more complex crystalline patterns of the LOFF condensate. Apparently the fast 
cooling of relatively young stars with a LOFF quark core is a consequence of the scaling laws for neutrino emissivity 
and specific heat. They depend on the existence of gapless points and follow from the existence of blocking regions in 
momentum space. Since this property is typical of the LOFF state, independently of detailed form of the condensate, 
a rapid cooling should be appropriate not only for the simple ansatz assumed in Eq. but, more generally, for any 
LOFF condensate. 
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APPENDIX A: QUASIPARTICLE DISPERSION LAWS AND MIXING COEFFICIENTS 



In this appendix we compute the quasiparticle dispersion laws, the gappless points and the Bogoliubov coefficients 
for the electrically and color neutral LOFF state of QCD with three flavors in beta equilibrium. In the basis A = 



(1,...,9) - K 



s r , Ub, s g , db) the gap matrix is as follows: 
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We consider in detail only the solution with Ai = 0, A 2 = A3 = A (notice that a factor e 2lq2 ' r = e 2lq3 ' r has been 
omitted in the entries A 2 , A3). The case of the LOFF2s model (Ai = A 2 = 0, A 3 ^ 0), can be treated in a similar 
way. Since the matrix (|A1() is block-diagonal, we consider its various sectors separately. 



1. Sector A = (u r , d g , s b ) 



To get the dispersion laws of the quasiparticles one has to find the poles of the propagator. This corresponds to 
solve the equation detS* -1 = 0. The procedure is simplified because one can write the determinant as a product of 
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various minor determinants. In particular for the quasiparticles given by linear combinations of u r , d g and Sf, quarks 
one has to evaluate the determinant of the matrix 



wi A A 
A ii 2 
A lj 3 

where 

wi = e — p + /j,i, lu 2 = e + p - ^2 , &3 = e+p- (J.3 . 
The dispersion laws are obtained by solving the cubic equation in the variable e 

(e-p + /ii)(e+p-/i2)(e+p-M3)-A 2 [2(e + p)- ( u 2 - M3 ] = 0. 
If we define yj = p — [ij (j = 1, 2, 3) and 

d = 2/12/22/3 + A 2 (?/2+ 2/3) 

c 



(A2) 



(A3) 



(A4) 



b = 2/1 - 2/2 - 2/3 
b 2 

P = C~ T , 



c = -ym - ym + 2/22/3 - 2A 

, be 2b 3 
q = a — - 



then the roots of Eq. I|A4|I are as follows 



eo 



3w 



pw 

T 



27 



C w 



27q 



VlV + 81g 2 



P 

C, w 



(A5) 



(A6) 



with w = 



One can show that only the quasiparticle of energy eo(p) is gapless, while e aj & refer to gapped 

quasiparticles (this is true for the reference value we use, y = 140 MeV, and for any y < 148 MeV; for y > 148 we 
find that are still gapped except for a tiny region around the North and the South poles). There can be one 
gapless mode (Pj) or three (Pi, Pn, Pm) depending on the angle that 2q forms with difference of the pairing quark 
momenta. The gapless modes are solutions of the equation 



(P - M1 )(P - M2)(P - /x 3 ) + 2A 2 (P 
The Bogoliubov coefficients given by 



M2 + M3 







B 2 Ur (0 - 



(Q 2 U>3) 2 

f(uj 2 ,0J 3 , A) 
(A^ 3 ) 2 

f(w 2 ,Ul3, A) 

(AQ2) 2 

/ (0)2,0)3, A) 



(A7) 

(A8) 
(A9) 
(A10) 



with f(ai, ct2, CK3) = af + afa§ + a|a| , represent the probability that the gapless quasiparticle is respectively u r , 
d g or S(,. We note that the dispersion laws are anisotropic as they depend on the angle $ that the quark momentum 
form with the pair momentum. In Fig. 0|we plot the dispersion law for the gapless quasiparticle for cos^ = 0. The 
quark velocities Vi, vn and vm are given by the slopes of the curve near the gapless points. 



2. Sector A — (d r ,u g ) 
In this sector quasiparticles have dispersion laws as follows: 



I (d r ,U g 
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FIG. 4: (Color online) Dispersion laws for the quasiparticles given by linear combinations of u r , d g and s& quarks. One of the 
three laws is gapless (red online) with three gapless momenta. The graph is obtained at p = 500 MeV, m,//! = 140 MeV and 
for cosf? = 0, where cosi? is the angle between the quark momentum and the pair momentum 2q. 
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FIG. 5: (Color online) Solid and double dotted-dashed lines denote gapless dispersion laws (red and black online, respectively 
for u g — d r and Ub — s r ), dashed and dot-dashed lines refer to gapped dispersion laws (green and blue online, respectively for 
u g — d r and Ub — s r ). Values of the parameters as in Fig. 0] 



where z = cost? and i? is the angle that the quark momentum form with the pair momentum. Since 

, in - Ms y y {, z\ , A10 s 

d ^^^-8- qZ ^8{ 1 ~-J ' (A12) 
with y given in Eq. Gapless modes are present for 

^ 2 = (| - <iz) 2 > A 2 , (A13) 

therefore the existence of gapless momenta depends on the value of z = cost?. The dispersion law e + is gapless for 
S/i < whereas e_ is gapless in the other case. The corresponding gapless momenta are 



,2 / - X 2 
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Denoting the gapless dispersion laws with e g , the mixing coefficients are given by 

= e g -p + K = A 

Their values at gapless momenta are 

Bl(P^) = Bl g {P u + ^) = \(l- VS ^ A2 ^j , (A16) 

Bl(P^) = Bl(P^) = \ [x + . (A17) 

We notice that 8^ > corresponds to z < z q a condition that is true almost everywhere (except in a tiny region 
around the North pole). 

3. Sector A — (s r , Ub) 
In this sector the quasiparticles have dispersion laws 

I J.S r ;U b ) 



In this case 



Gapless modes are present for 



The gapless momenta are 



'-q 

The values of the Bogoliubov coefficients at the gapless points are 



Now S/i > corresponds to z < ~z q . 



4. Sector A = (s g , db) 



There is no mixing in this sector because Ai = 0, therefore quasiparticles have dispersion laws 



e 



Gapless momenta are 



(A18) 



fj, 6 -Hr y y z . 

^=—2— = "8 ~ qZ = -*{ 1+ *) ■ (M9) 



b[? = (| + > A 2 . (A20) 



^^ ± Jy-( 1+ ±) , (A2 i) 



B 2 A pu bSr) = B l b{ pu b s r) = 1 L_ y/W_ AA (A22) 



Bl(P?*) = Bl b (P^) = \ h+ VS ^ A j • (A23) 



\P-V8,9\- (A24) 



P Sa = M8 = M-§, (A25) 
P d > = M9 = M+^ • (A26) 
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The mixing coefficients are B Sg (P Sg ) = Bd b {P db ) = 1 and the corresponding velocities are equal to unity. 
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